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AePW building block approach to validation 


Utilizing the classical considerations in 
aeroelasticity 

• Fluid dynamics 

• Structural dynamics 

• Fluid/structure coupling 



Unsteady fluid dynamics 



AePW-1: Focused on Unsteady fluid dynamics 


AePW-2: Extend focus to coupled aeroelastic simulations 






You are invited to participate in AePW-2 


Extend focus to coupled aeroelastic simulations 
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AePW-1: Applying the Lessons Learned 


• One configuration only 

• Benchmarking case: including a case that we have confidence can be 
“well-predicted” 

• Comparison metrics: 

- Unsteady quantities for all cases 

- Integrated sectional forces and moments 

- Critical damping ratios and frequencies 

- Extended statistics: mean, std, mode, max, min 

• Time histories from solutions requested because 

- nothing is steady 

- single person, single method of post-processing matters 

- there’s always more to see- nonlinearities, off-nominal frequency content 

• Results requested at more finely spaced points than experimental data 

• Common grids suggested for analyses 

• Various fidelity aerodynamic contributions encouraged 

• Discussion telecons for analysis teams 



Overview of requested submittal data sets 


• Steady rigid pressure coefficient distributions: 
statistics of the results 

• Time histories 

o Angle of attack 

o Leading and trailing edge displacements 
o Pressure coefficients 
o Lift & pitching moment coefficients 
o Sectional lift & pitching moment coefficients 

• Frequency response functions: C p /0 

o At forced oscillation or flutter frequency 
o Across 0-100 Hz 

• Static aeroelastic pressure coefficient 
distributions: statistics of the results 

• Flutter bounds 



Example simulation 
results: 

Slices through the pressure field at 
different points in the flutter cycle 

#1 





Cp vs time and x/c 
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The results shown here are at 
60% span for the upper surface 
at the dynamic pressure of the 
experimental data set (168.8 
psf) 

Data sets are requested at 

• 60% and 95% span 

• Upper & lower surfaces 

• Experimental flutter 
condition (168.8 psf) & 
Computational flutter 
condition as determined 
by each analysis team 






Frequency, Hz 


Example Flutter 
summary results 
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These example results were 
calculated using 
o URANS + SA 
o Medium fidelity grid 
o Relatively coarse time step 
The generalized displacement time 
histories were analyzed to produce 
the damping and frequency results 



Steady rigid 
pressure 
distributions 


For the primary forced oscillation case, Case #1 , disagreements with experimental 
data limited to the peak of the upper surface shock. 

For the primary flutter case, Case #2, shows a well-matched rigid pressure 
distribution without much variation among the computational results. 

The complexity of the Case #3 is indicated by the variation among the 
computational results & difference from the experimental data ->Shock location, 
shock strength, aft loading especially on lower surface. 


Case comparisons 




Temporal parameter influences 
on aeroelastic stability results 


FUN3D analysis (URANS + SA) 




Case 1 


Optional Case 3 
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Summary of temporal parameters 
for different solutions 



Summary of temporal properties 
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Varying time step size at q = 168.8 psf; 25 
subiterations per global time step 


BSCW FUN3D URANS + SA Dynamic Aeroelastic Analysis 
Medium Grid 

Mach 0.74, Mean angle of attack 0 degs, 

Dynamic pressure 168.8 psf 
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Varying time step size. Medium Grid 
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Mode 1 
Mode 2 
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DT = 1.2, dt = 0.0002 sec sample (Fine) 

DT = 12.1, dt = 0.002 sec/sample (Medium) 
DT = 24, dt = 0.004 sec'sample (Mod2) 

DT = 61, dt = 0.01 sec'sample (Modi) 

DT 121, dt = 0.02 sec/sample (Coarse) 


Mode 2 
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For constant number of subiterations, 25: 

• Simulation shows stable results for 

DT >= 24 (refinement later showed 
That the neutrally stable DT at this 
dynamic pressure is 21 .2?) 

• Smaller time steps show unstable behavior. 

• Larger time steps show stable behavior. 





Critical damping ratio. 


Varying time step size at q = 168.8 psf; 

Temporal error convergence 10%, 

1000 subiterations maximum per global time step 


Stability at q=169 psf, Mach 0.74, a= 0° 
Varying time step size. Medium Grid, TC 10% 
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Physical time step size, seconds 

For temporal error convergence of 1 0%, with 
gvelO = 5.0 on both modes: 

• Simulation shows unstable system for all cases 
gvelO = 0.5 on both modes: 

• Simulation shows stable Mode 1 behavior at smallest time step 
(more iterations running to see if this changes) 


Jen, remember that you are 
Assuming that there are 2 eigenvalues, 
One stable and one unstable. 

The fine grid result with gvel = 0.5 may just 
Be indicating this other root. 

Need to combine the gvel results and then 
analyze the resulting signal for 2 modes. 


Smaller time step = more unstable 
Higher value of initial kick = more unstable 


Critical damping ratio. 


Varying time step size at q = 168.8 psf; 

Temporal error convergence 10%, 

1000 subiterations maximum per global time step 

Stability at q=169 psf, Mach 0.74, a= 0° 

Varying time step size. Medium Grid, TC 10% 

ATheta, gvel 5.0 
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♦ Mode 1 

♦ Mode 1, gvel 0.5 
□ Mode 2 
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For temporal error convergence of 1 0%, with 
gvelO = 5.0 on both modes: 

• Simulation shows unstable system for all cases 
gvelO = 0.5 on both modes: 

• Simulation shows stable Mode 1 behavior at smallest time step 
(more iterations running to see if this changes) 

• Smaller time step = more unstable 

• Higher value of initial kick = more unstable 


SAME INFORMATION 
AS ON PREVIOUS 
SLIDE, but showing 
Horizontal axis as 
Time steps/4 Hz cycle 


The following 4 slides are from 
Pawel and show the results for 
Coarse, Medium and Fine grids for 
DT = 1, with 10% temporal error 
convergence and 1444 as the 
maximum number of subiterations 



Current FUN3D results: 
spatial and temporal convergence 
AePW-2 Case#2 Flutter results 
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Coarse Grid, q = 168.8 lb/ft A 2 
DT = 1, 10% Temp. Error 
1444 timesteps per cycle 
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Current FUN3D results: 
spatial and temporal convergence 
AePW-2 Case#2 Flutter results 
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Medium Grid, q = 168.8 lb/ft A 2 
DT = 1 , 10% Temp. Error 
1444 timesteps per cycle 
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Current FUN3D results: 
spatial and temporal convergence 
AePW-2 Case#2 Flutter results 




Current FUN3D results: 
spatial and temporal convergence 
AePW-2 Case#2 Flutter results 


DT = 1, 1444 time steps per cycle, 

10% Temporal error controller, q = 168.8 lbf/ft A 2 




Varying time step size at q = 168.8 psf; 

Temporal error convergence 10%, 

1000 subiterations maximum per global time step 

Damping values were calculated using the generalized displacements associated 
With the two aeroelastic modes, treating each as if they contained only 
A single mode. Near neutral stability, this isn’t a bad assumption. 

How good are the fits that are used in the damping calculations? 

The following plots show this for the Mode 2 data sets. 



gdisp 


Q 169, DT 1.2, comparison of temporal error 
convergence & fixed number (25) subiterations 



q = 168.8, DT = 1.2, Temporal Error Convergence 
Criteria 10%, Maximum # of subiterations: 1000 



Q 168.8 psf, DT 24 comparison with DT=1 10% 

temporal error convergence 


Dl 24, Temp Error 0.02, mode 2 
Dt 24, Temp Error 0.05, mode 2 
Dt 24, Temp Error 0.1 0, mode 2 
Dt 24, Temp Error 0.20, mode 2 
DT 24, 25 subiters, mode 2 
DTI. 2, Temp Error 0.10, mode 2 


BSCW analysis in FUN3d + SA 
Medium Grid 

Mach 0.74, qbar 1 69, aoa 0 
Flutter analysis 
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Changing the subiteration criteria after stability 
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BSCW analysis in FUN3d + SA 
Medium Grid 

Mach 0.74, qbar 169, aoa 0 

Flutter analysis with time step DTmod2, 24.375 


behavior is established 



I I L 

650000 700000 750000 

Simulation Time 


BSCW analysis in FUN3d + SA 
Medium Grid 

Mach 0.74, qbar 169, aoa 0 

Flutter analysis with time step DTmod2, 24.375 
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gdisp 


Varying temporal error convergence criteria 


Dt 24, Temp Error 0.02, mode 2 
Dt 24, Temp Error 0.05, mode 2 
Dt 24, Temp Error 0.10, mode 2 
Dt 24, Temp Error 0.20, mode 2 


BSCW analysis in FUN3d + SA 
Medium Grid 

Mach 0.74, qbar 1 69, aoa 0 
Flutter analysis 



The damping increases 
as the temporal error criteria 
becomes less stringent. 

(i.e £ increases with increasing 
temporal error convergence 
percentage) 

System is stable for 20% TEC 

All other solutions show flutter 
(instability) and then limit cycle 
behavior. 

Limit cycle onset occurs at different 
simulation time for each case. 

I have the subhist and fun3d.out 
files for the 2% and 5% TC cases. 
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Flutter solution starting from rigid solution vs 

static aeroelastic solution 


(DT = 24.375, 5% temporal error convergence) 


BSCW analysis in FUN3d + SA 
Medium Grid 

Mach 0.74, qbar 169, aoa 0 

Flutter analysis with time step DTmod2, 24.375 
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Mode 1, initiated from static aeroelastic solution 
Mode 2, initiated from static aeroelastic solution 
Mode 1, initialed from rigid solution 
Mode 2, initiated from rigid solution 
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Solution initiated from 
The rigid solution shows 
table behavior. 

Solution from static 
aeroelastic solution 
shows unstable 
behavior and then 
limit cycle oscillation 


Physical LCO prediction? 
Q 169, DT 15, TC 10% 


BSCW FUN3D Flutter Analysis, Medium Grid 

Mach 0.74, 0 degs AOA, qbar 168.8 psf 

DT = 15, dt = .0025 sec/sample, 100 samples/4Hz cycle 

Temporal Error Convergence 10%, 1000 maximum subiterations 
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Initial velocity kick (gvelO) variations 


For case 1250 time steps/cycle (DT = 1.2), 
Medium grid, 10% temporal error convergence, 
qbar = 168.8 psf 

• Currently running 0.5, 2.75 & 10.0 


For case ~ 200 times steps/cycle (DT = 7), 
Medium grid, 10% temporal error convergence, 
qbar = 168.8 psf 

• Ran gvelO = 5.0 & gvelO = 0.5 



Velocity kick influence 
Q = 168.8, DT = 7 
TC 10% 



Velocity kick influence 
Q = 168.8, DT = 1.2, TC 10% 
Partial results. Dec 8, 2014 

BSCW Flutter Analysis using FUN3D (URANS + SA), Medium Grid 

qbar 169 psf, DT =1.2, 10% Temporal Error Convergence, 1000 maximum subiterations 



Dec 8, 2014 

Thoughts regarding LCO results 

• If these are physical LCO results, then 
regardless of the gvel, the results should go to 
the same magnitude? 

• That is, if they do not encounter some 
violation or explosion due to numeric 

• Hmmm. Should they? Or, if it's physical, 
shouldn't the size of the velocity perturbation 
influence the results? Basins of attraction and 
all that? 



gdisp 


152 psf. 

Varying DT and subiteration convergence specification 


BSCW Flutter Analysis, Medium Grid 

Mach 0.74, 0 degs aoa, qbar 1 52 psf dt=i2, tc o.io, mode i 

DT=24,TC 0.10, model 

DT=1.2, 15 subiter, mode 1 
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BSCW Flutter Analysis, Medium Grid 
Mach 0.74, 0 degs aoa, qbar 152 psf 
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135 psf, 

Varying DT and subiteration convergence specification 


BSCW Flutter Analysis, Medium Grid 
Mach 0.74, 0 degs aoa, qbar 135 psf 


DT=12, 25 subiter, mode 1 
DT=12, TC 0.05, model 
DT=1.2, 25 subiter, mode 1 
DT=1.2, TC 0.10, model 
DT=1.2, TC 0.02, mode 1 
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BSCW Flutter Analysis, Medium Grid 
Mach 0.74, 0 degs aoa, qbar 135 psf 
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DT=12, 25 subiter, mode2 
DT=12, TC 0.05, mode2 
DT=1.2, 25 subiter, mode2 
DT=1.2, TC 0.10, mode2 
DT=1.2, TC 0.02, mode2 
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